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Abstract 

Luminosities and fluxes of the expected population of galactic gamma-ray 
pulsars become foreseeable if physical distributions at birth and evolutive his¬ 
tory are assigned. In this work we estimate the contribution of pulsar fluxes 
to the gamma-ray background, which has been measured by the EGRET 
experiment on board of the CGRO. For pulsar luminosities we select some 
of the most important gamma-ray emission models, taking into account both 
polar cap and outer gap scenarios. We And that this contribution strongly 
depends upon controversial neutron star birth properties. A comparison be¬ 
tween our simulation results and EGRET data is presented for each model, 
finding an average contribution of about 10%. In addition, we perform the 
calculation of the number of new gamma-ray pulsars detectable by GLAST 
and AGILE, showing a remarkable difference between the two classes of mod¬ 
els. Finally, we suggest some improvements in the numerical code, including 
more sophisticated galactic m odels and different populations of pulsars like 
binaries, milliseconds, anomalous pulsars and magnetars. 


1 Introduction 

The estimation of the contribution of pulsars to the gamma-ray background is 
based on the comparison between theoretical predictions on the flux emitted by 
these sources and experimental observations. Each model for the emission of 
radiation by a stellar population is specified by the galactic distribution of sources 
and their individual luminosities. 
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Figure 1. Distributions of pulsars ages (t), distances (d), magnetic fields (Bq), periods (Po)i 
period derivatives (P). Dot lines represent a situation with constant magnetic field, solid lines 
are referred to the field decay case (E. Bisesi, 2002). 


2 Galactic model 

2.1 SPATIAL DISTRIBUTIONS 

We here present a numerical simulation of the galactic population of radio pulsars 
giving all initial spatial and physical parameters in the framework of the galactic 
models of Paczyhski (1990) and of Gonthier & al. (2002) we then make them 
evolve accordingly with the velocity model of Stumer & Dermer (1994) ^H]- We 
show histograms of distributions of various pulsar properties in Fig. o, assuming 
both constant magnetic field and field decay. 

2.2 PULSAR LUMINOSITIES 

In the commonly accepted hypotesis that every pulsar irradiates on the whole 
electromagnetic spectrum, we assign a gamma luminosity to each source. We 
consider a suitable choise of emission models. 

POLAR CAP MODELS: 

• Harding (1981) gl: L.^(>100 MeV) = 1.2 x ph s"!; 

• Zhang & Harding (2000) [1 2| : 

L.^(/) = 5.87x1035Ri 2®^^P"^/^ ph s-i for > 6.0, 

L^ill) = 1.0 X 1035Ri 2 ph s-i for p-ii/28 < 5 Q; 
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Figure 2. Map of the gamma-ray sky as seen by EGRET (E. Bisesi, 2002). 
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Figure 3. Comparison between the EGRET gamma-ray background (in blue) and predictions for 
the polar cap model of Harding (1981) (in violet). We also give the simulated number of pulsars 
above the EGRET detection threshold for this model. Plots on the right side are enlargements 
of the simulated profiles (E. Bisesi, 2002). 


• stumer & Dermer (1994) |10|: = 6.25 x P ^ ph s 


OUTER GAP MODELS with the death line 5 logB — 12 logP < 72: 

• Romani &: Yadigaroglu (1995) 12]: = 1.56x phs“^; 

• Cheng & Zhang (1996) ^ : = 3.93 x 10^"^Bi 2 °-^ p-°-^ ph s-h 

We calculate the integrated flux on the galactic latitude —10° < b < 10° as 
a function of the longitude 1; results are plotted together with the experimental 
points as they have been measured by the EGRET experiment. 

Fig. @ and respectively show the map of the gamma-ray sky as seen by 
EGRET and our results for the polar cap model of Harding (1981). 

The contribution of pulsars to the gamma-ray background is less than 10% for 
this model; the foreseen number of sources above the EGRET detection threshold 
of 10“^ ph cm“^ s“^ in not far from the observed value. 
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Figure 4- Diagram period—magnetic field for different pulsar families. Points indicate young 
isolated pulsars, circle—dots pulsars in binary systems. The population of millisecond pulsars is 
visible at the bottom left (from Phinney & Kulkarni, 1994 |^). 

3 Comparison among models 

In the previous discussion we have forced the physical parameters involved in the 
calculation of gamma ray luminositied and fluxes to assume some definite values. 
In order to give a meaningful comparison among different emission models an 
extension to a more general parameter space is required. 

So we define a parameter space whose extremes are corresponding to the limit 
conditions for pulsating gamma neutron stars; we then follow the evolution of a 
representative point inside this space. The most relevant physical parameters are 
the rotation period, the magnetic field and the velocity fron the Galactic Centre. 

3.1 PULSAR POPULATIONS 

Different families of pulsars follow distinct evolutive histories, and this has a di¬ 
rect impact on their physical properties. Furthemore, physical mechanisms ruling 
over the gamma-ray emission are peculiar of each class. Fig. 0 shows a dia¬ 
gram period—magnetic field for young isolated, binaries and millisecond pulsars. 
Magnetars are characterized by very high magnetic fields, typically 10^^ — 10^® G. 
We bind the parametric space in this way: 

• Pomin = 0.01 s (when the centrifugal force equals molecular bonds of the 
star, yielding its survival); 

• Pomax = 1 s (fixed from observations); 

• Bomin = 10^^ G (the lowermost value for young isolated pulsars); 

• Bomax = 4.413 X 10^^ G (condition for the magnetar regime); 
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Figure 5. Comparison among models (E. Bisesi, 2002). 


• vomin = 120 km/s (Gonthier & al., 2002 ^); 

• vomax = 350 km/s (Lyne & Lorimer, 1994). 

We show a tridimensional representation of the parametric space (logPo); 
(logBo), V in the left side of Fig. (Ol. The bidimensional portion on the right 
shows the positions of two models of interest nini. 

We compare predictions for the five models introduced above for a set of points 
of the parametric space. As physical mechanisms for the gamma-ray emission differ 
from one pulsar family to another, we would give luminosity function expressions 
for each of them. We purpose to improve our work in such a way in future; in this 
context we restrict our analysis to young isolated pulsars. An example relevant to 
the point is shown in Fig. ®. 
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Figure 6. Comparison between EGRET observations and simulation results for the point Xz of 
the parametric space (E. Bisesi, 2002). 



























6 


Contribution of pulsars to the gamma-ray background 


4 Discussion 

For each point considered, we evaluate the percentual excess between the foreseen 
number of gamma-ray pulsars for each model and their actual number of 7. The 
total percentual excess for each model gives us the possibility of select the most 
reliable model, which we find being that of Harding (1981) [H]. 

We finally estimate the number of gamma-ray pulsars detectable by the next 
space telescopes GLAST and AGILE for the model selected. 

GLAST and AGILE detection thresholds are 6 x 10“® ph cm“^ s“^ and 10“^ 
ph cm“^ s“^ respectively. 

Our results for the five points of the right side of Fig. 0 are shown in Table I. 

Predictions for GLAST are very optimistic, remarkably we expect to detect a 
very large number of new gamma-ray pulsars, opening very promising frontiers in 
understanding these mysterious and fascinating objects. 


Table 1. Total percentual excess for each model on the whole parametric space for young pulsars. 


Model [Total percentual excess] = [Percentual excess(X^)| 


H 

27 

ZH 

220 

SD 

174 

RY 

30506 

CZ 

2277 
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